Modelling one-dimensional driven diffusive systems by the Zero-Range Process 
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The recently introduced correspondence between one-dimensional two-species driven models and 
the Zero- Range Process is extended to study the case where the densities of the two species need not 
be equal. The correspondence is formulated through the length dependence of the current emitted 
from a particle domain. A direct numerical method for evaluating this current is introduced, and 
used to test the assumptions underlying this approach. In addition, a model for isolated domain 
dynamics is introduced, which provides a simple way to calculate the current also for the non-equal 
density case. This approach is demonstrated and applied to a particular two-species model, where 
a phase separation transition line is calculated. 



I. INTRODUCTION 

Phase separation in one-dimensional driven systems 
has attracted much attention of late Q, 0, B|- In 
contrast to equilibrium one-dimensional systems, where 
phase separation cannot occur unless the interactions 
are long ranged, several examples of phase transitions 
in one-dimensional non-equilibrium steady states have 
been given 

HE ISSUES EH- These models generally 
have local noisy dynamics and some conserved quantity 
or quantities driven through the system. 

Particular attention has been paid to a simple yet gen- 
eral class of models with two species of particles which 
are conserved under the dynamics P, Q . These models 
are defined on a ring, where each site can take one of 
three states: vacant, occupied by a positive particle, or 
occupied by a negative particle. Two conservation laws 
are obeyed by the dynamics, which can be cast into two 
conserved quantities - the total density of particles in the 
system, p, and the fraction r\ of positive particles out of 
the total number of particles. 

To study these models a coarse-grained description has 
been developed [j| . In this description one views the mi- 
croscopic configuration of the model as a sequence of par- 
ticle domains, bounded by vacancies. Each domain is de- 
fined as a stretch of particles of both types. The idea is to 
view particle domains as urns which may exchange parti- 
cles. At a coarse-grained level one identifies the current of 
particles through a domain as the hopping rate of parti- 
cles between neighbouring urns (see Fig.^). This coarse- 
grained description then defines a Zero-Range Process 
(ZRP) for which the steady state may be solved exactly. 
When such correspondence is applicable one can use the 
ZRP to obtain the distribution of the domain size, al- 
though information about the correlation between the 
two species of particles is lost. 

Generally the identification of the driven system with 
a ZRP is at a coarse-grained level and is not exact. 
Rather, it relies on the applicability of some physical as- 
sumptions, as discussed below. However, for a particular 
model for which an exact solution of the steady state ex- 
ists, it could be shown that the mapping of the steady 
state to that of a ZRP is indeed exact . 
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FIG. 1: A microscopic configuration of a two-species driven 
model (bottom) and its corresponding configuration in the 
ZRP (top). 

The correspondence between driven models and the 
ZRP may be used to address the question of existence of 
phase separation in the driven model . In such a tran- 
sition a fluid phase, where the distribution of particles 
and vacancies is homogeneous, becomes upon increas- 
ing p a phase separated state. This state is character- 
ized by macroscopic particle domain devoid of vacancies. 
In the context of the ZRP the transition into the phase 
separated state corresponds to a condensation transition 
whereby one urn becomes filled by a macroscopic number 
of particles HG3- 

A careful analysis of the corresponding ZRP suggested 
the following criterion for phase separation. Let Uk be 
the rate at which particles leave a domain of size k. If 
Uk vanishes in the thermodynamic limit k — > oo, phase 
separation is expected at any density. Otherwise, for 
large domains, Uk typically takes the form 

Uk = "oo ^1 + jrj , (1) 

where and b are constants. The existence of phase 
separation is then related to the value of b. Phase sepa- 
ration cannot exist as long as b < 2. However, if b > 2 
a phase transition into a phase separated state occurs as 
the particle density is increased. This picture has been 
used to argue against the existence of phase separation in 
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some models [9j , to study crossover phenomena |13J , and 
to suggest a model which does exhibit a phase transition 
into a phase separated state [Icj . 

In order to apply the ZRP picture for a given model 
one has to evaluate the rates Uk by which domains ex- 
change particles. This may be too difficult a task to carry 
out analytically, and may involve exceedingly large com- 
putation time to estimate numerically. However, it has 
been suggested that in order to estimate Uk one may re- 
duce the full many-domain system into a single isolated 
domain problem. This is done by modelling an isolated 
domain by an open chain of particles exchanging particles 
with reservoirs at its ends. At the boundaries, particles 
are injected and ejected at constant rates. When these 
rates are large enough, the system is known to be forced 
into a maximal-current phase with a bulk density 77 = 1/2 
[l4j . This approach is therefore only applicable when the 
densities of the two species of particles in the full model 
are equal, and each domain is stationary on average. In 
this case the rate Uk is directly related to the current Jk 
which flows through the domain. Thus, for equal densi- 
ties of positive and negative particles the many-domain 
problem is simplified to the problem of a single domain 
of fixed size. The calculation of the current through the 
single domain ma y b e tackled numerically or analytically, 
when possible |9l- flO| - 

It is important to notice that the approach outlined 
above for describing a many domain system by a ZRP 
relies on two assumptions. Firstly, it is assumed that 
domains are uncorrelated, in the sense that the current 
Jk flowing through a given domain depends only on its 
own size. Secondly, the current through a domain of 
length k is assumed to take its steady state value with 
respect to a system of size k (even though k fluctuates) 
and this steady state value is identified with that of an 
open system. 

In this paper we seek to test these underlying assump- 
tions. Whereas in previous studies the current Jk of a 
domain of length k was studied using an isolated single 
domain, here we introduce a numerical method for di- 
rectly measuring Jk of a fluctuating domain within the 
full system. An agreement between the two methods val- 
idates the assumptions behind the criterion. 

We also seek to extend the approach to treat the case 
of non-equal densities. As discussed above, modelling a 
single domain by an open system is not applicable in this 
case. Moreover since we have unequal densities of posi- 
tive and negative particles, the vacancies are not station- 
ary and drift on average. Thus in order to reduce the 
many-domain problem to a single domain problem we 
have to deduce the appropriate ensemble for the single 
domain. In this paper we propose a model for a fluctuat- 
ing isolated single domain which can be used to calculate 
a steady state current Jk in a domain of non-equal densi- 
ties. This result is checked by using the direct measure- 
ment discussed in the previous paragraph. Moreoever in 
a limiting case we can solve the model exactly through 
a matrix product ansatz and show that it produces the 



correct ensemble. 

This paper is organized as follows. In Section 2 we 
define a specific model that we use to demonstrate our 
approach. In Section 3 we describe the method for direct 
measurement of the currents within the many-domain 
system, and apply it to both cases of equal and non- 
equal densities. A model for an isolated domain which is 
not restricted to the case of equal densities is introduced 
in Section 4 and an exact solution of a limiting case is 
given. We then discuss generalizations of the approach 
to other models in Section 5, and present a summary and 
outlook in Section 6. 



II. MODEL DEFINITION AND SOME KNOWN 
RESULTS 

In order to study the correspondence between driven 
diffusive systems and the ZRP, we consider in the main 
body of this paper a particular driven model as a test 
case. In this section we define the model and present 
some analytical results. 

The model is defined on a one-dimensional ring of L 
sites. Each site i is associated with a 'spin' variable Si. A 
site can either be vacant (si = 0) or occupied by a posi- 
tive (si = +1) or a negative (si — —1) particle. Particles 
are subject to hard-core repulsion and a nearest-neighbor 
'ferromagnetic' interaction, defined by the potential 

V = -j£>*+i ■ (2) 

i 

Here < e < 1 is the interaction strength, and the sum- 
mation runs over all lattice sites. The model evolves ac- 
cording to the nearest-neighbour exchange rates 

+ - 1 ±^ y -+ +0^0+ 0-^-0, (3) 

where AV is the difference in the potential V between 
the initial and final states. The number of particles 
of each species, N + and AT_, are conserved by the dy- 
namics. Alternatively, the system can be character- 
ized by two conserved densities, namely the total den- 
sity p — (N + + N-)/L, and the relative density 77 = 
A^-i- / (A^ + N-) . This model, which is a generalization 
of the Katz-Lebowitz-Spohn model 0, , was intro- 
duced in and studied for the case of equal densities 
of positive and negative particles, rj = 1/2. 

A. The non-interacting case, e = 

Let us first discuss the case e = 0, where particles only 
interact through the hard-core exclusion. In this case 
an exact solution shows that within a grand-canonical 
ensemble to be defined below domains are uncorrelated, 
and the steady-state weight factorizes into a product of 
single-domain terms. This is the case for both equal and 
non-equal densities of the two species. These results are 
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obtained by considering a grand-canonical ensemble in 
which the number of vacancies M is kept constant. The 
number of particles of each species, and thus the size of 
the lattice, are allowed to fluctuate. A fugacity £ is at- 
tached to the positive particles, thus controlling the rel- 
ative density between the two species (£ = 1 corresponds 
to equal densities). All configurations of the system can 
be described in terms of domains of particles, where a do- 
main is defined as an uninterrupted sequence of particles 
of both species. The weight Wm ({ki}) of all configura- 
tions in which fej particles reside in the ith domain is then 
given by 



M 



w M ({h}) = l[z k >z ki (0 



(4) 



i=l 



where Zk is the sum over all weights of microscopic con- 
figurations of a domain of length fc, and z is the fugac- 
ity, which controls the overall density of particles in the 
system. Hence, within the grand canonical ensemble do- 
mains are statistically independent with a domain size 
distribution 



P(k) ~ z k Z k {£,) . 
The grand canonical partition function is given by 



A I 



(5) 



(6) 



{fc,}fc 



In the equal density case the exact solution reveals 
that Zk is identical to the partition function of the to- 
tally asymmetric exclusion process on a one- dimensional 
lattice of k sites with open boundary conditions . For 
non-equal densities, it turns out that Zk is identical to 
the grand-canonical partition function of the same pro- 
cess defined on a ring of size h + 1 with a single vacancy 
[H Ell- In both cases 
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for large k. The resulting distribution function (JjjJ im- 
plies that the model does not exhibit phase separation 
at any density. In particular, for any choice of p one can 
choose the fugacity z < 1 such that the average density 
satisfies p/(l — p) = J kP(k)dk. 

Moreover, in the case e = the correspondence be- 
tween the steady-state of the model and that of the ZRP 
can be made explicit. Within a grand-canonical ensemble 
of the ZRP urns are statistically independent. The dis- 
tribution function P ZRP (fc) for the occupation of a single 
urn is given by P ZRP (fc) ~ z k Yl m — 1 l/u m . On the other 
hand, in the e = driven model the steady-state current 
Jk flowing through a domain is given by Jk = Zk-i/Zk- 
The steady-state distribution function (JjjJ can then be 

written as P(k) ~ z k l\t=i 1 l J m- Thus P(k) = P ZRP (fc) 
with Uk = Jk- Using Q one obtains 6 = 3/2 for this 
case, implying no phase separation. 



B. The interacting case, e > 0, at 77 = 1/2 

We now turn to the more general case, e / 0. Here 
no exact mapping to the ZRP is available. However, it 
was conjectured in 0, 0] that the physical picture ob- 
tained for the non-interacting case remains valid, namely 
that the hopping rates of a corresponding ZRP should 
be identified as the steady-state currents of isolated do- 
mains. One therefore needs to calculate the asymptotic 
form of the steady-state current running through a do- 
main, Jfc(e) ~ Joo(e) (1 + 6(e)/fc). For r) = 1/2, where the 
average domain velocity vanishes, the current Jk may be 
calculated by considering an isolated domain with open 
boundaries, which exchanges particles with reservoirs at 
its ends at high rates. 

It has been argued [20I l2l| that the coefficient 6(e) of 
an isolated open domain is given by 



b(e)=cb R (e,r) = 1/2) 



(8) 



Here 6#(e, 77) is the coefficient corresponding to a closed 
fully-occupied ring, and c is a universal constant which 
is equal to 3/2. For a ring the coefficient bjt(e, rj) can be 
calculated at any density 77 [ljj. It is given by 



A(e, 77)^(6,77) 
2Joo(e, rj) 



(9) 



In this expression A = d 2 J^/drj 2 . The compressibility 
k = linifc^oo k^ 1 [(rj 2 ) — (rj) 2 ) is evaluated in a grand- 
canonical ensemble of a fully-occupied ring with average 
density (77) = 1/2. Using the known properties of the 
steady state of this model [la . Hq it can be shown that 



(10) 



K(e,rj) = 77(1-77)^/1 + 477(1-77) 



1 + e 
1 -e 



1 11) 



(7) where 7 = [477 (1 - 77)]- 1 / 2 + [(4?7(l-?7))- 1 -l+(l+e)/(l- 



^l 2 . Inserting (fTUjl and lfTT|) into JjJ) with r) = 1/2 one 
obtains the coefficient 6 of an open domain as a function 
of e for the equal density case. It is found that 6 > 2 
for e > 4/5. The criterion mentioned above therefore 
implies that in this model phase separation takes place 
at high density p for any e > 4/5. Note that although 
the expression of 6^(6,77) is valid for arbitrary 77 in the 
ring geometry, the resulting 6(e) (Eq. [HJ is relevant to a 
domain with open boundaries only at 77 = 1/2. A single 
domain model which is applicable for 7/ ^ 1/2 will be 
discussed in Section llVl 



III. 



DIRECT NUMERICAL MEASUREMENT 
OF DOMAIN DYNAMICS 



We now describe our method for a direct numerical 
measurement of the flow of particles Uk out of a domain 
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FIG. 2: Density profile of a domain of size k = 100 at equal 
densities (77 = 1/2) and e = 1/2. Solid line corresponds to 
direct simulation of the full many-domain model. Plus signs 
correspond to an open isolated domain, and squares corre- 
spond to the non- conserving single-vacancy ensemble. 
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FIG. 3: Results of numerical simulations of the full model for 
■q — 1/2 and p = 0.7, from which 6(e) and Moo are estimated. 
Solid lines are obtained by linear fits of the data. 



of size k. As discussed in the introduction, this method 
involves the full many domain system, and thus allows a 
check on the validity of reducing to the problem to that 
of a single domain. 

The idea is during a simulation of duration T to record 
the number of hopping events out of a domain of size k, 
and the average number of domains of size k. The ratio 
of these quantities yields Uk- Here one unit of time cor- 
responds to a single Monte-Carlo sweep. We thus define 



to+T 



to+T 



N, 



i £ m k (t) F k = ±J2 AW 



(12) 



t=t 



Here mk is the number of domains of size k residing in the 
system at time t, and fk(t) is the number of exchanges 
+0 — * 0+ occurred between times t and t + 1 at the 
boundaries of domains of size k. Of course, one can also 
define fk through the transition rates 0— — ► —0, with- 
out changing the following discussion. The measurement 
starts at time to, after short-time relaxations are over. 
Clearly, 



u k = Jim — . 

T->oo Nk 



(13) 



In practice, the measurement time T taken to be large 
enough to ensure convergence. Estimates for and b 
are then obtained from the linear fit of Uk to 1/fc. How- 
ever one can exploit the data obtained from numerical 
simulations better by integrating the distributions. Thus 
we define 



l=k 



i=k 



— ~-i ^-^ No 
Q k=Nk Y,~f 



i=k 



Using Uk = Uoo(l 4- b/k) one has 

U k = r=^r — Moo f 1 + bQ k 

N k V 



and one obtains b and Unn from a linear fit. 



(14) 



A. Equal densities, rj = 1/2 

While simulating the dynamics of the full model we 
have recorded the density profiles of domains of a give 
size k. These profiles are compared with those obtained 
from the single open domain calculation in Fig. [21 Both 
profiles are identical to within the statistical fluctuations 
except at sites 1 and k where there are systematic de- 
viations. The deviations at these two sites are to be 
expected as this is where the dynamics for the single do- 
main is simplified from the full model. The excellent 
agreement between the profiles indicates that the single 
open domain properly models a fluctuating domain in the 
full system. 

We now apply the direct numerical method outlined 
above to the case 77 = 1/2. The results of numerical 
simulations are given in Fig. for several values of e. 
Table[I]summarizes the resulting parameters 6(e) and Moo- 
For e = one knows from the exact correspondence of the 

model to ZRP that Uk — Jk = \ (l + ^ff^j ■ O ur direct 
measurement of v,k recovers these results quite faithfully. 
In general, for e > we find that the measured values 
of b and Moo are in close agreement with (IH Hllll . Thus 
we conclude that the current flowing through a domain in 
the model can indeed be viewed as the stationary current 
flowing trough an isolated open system of the same size. 
The proposition that by increasing e a phase transition 
into a phase separated state occurs is thus verified. 



B. Non-equal densities, r\ 7^ 1/2 

We now consider the model with non-equal densities 
of the two particle species, i] ^ 1/2. In order to apply 
the criterion for phase separation one needs to calculate 
the coefficient b(e,ri). As mentioned in Section [H] an 
exact mapping to the ZRP is only available at e = 0, 
where one finds b = 3/2 for any value of 77. To test 
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measured 





1.5 


1.51 


0.25 


0.25 


0.25 


1.67 


1.65 


0.2113 


0.2115 


0.5 


1.82 


1.79 


0.1585 


0.1583 


0.8 


2 


2.05 


0.075 


0.0749 



TABLE I: Results of direct measurement of Uk in the full 
model, compared with analytical results of Section 2. We 
estimate the error in the measurement of itoo to be 0.0005 
and of b to be 0.05. 
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FIG. 4: The coefficient & as a function of the density rj of 
positive particles within a domain. Solid lines are given by 
JSJ, crosses are direct measurements in the full system, while 
circles correspond to the non-conserving vacancy ensemble. 



the validity of the correspondence to ZRP and evaluate 
b(e, rj) for e ^ and arbitrary r\ one has to resort to direct 
numerical simulations, as discussed above. 

We have simulated the model for e = 0.25, 0.5 and 
r] = 0.6, 0.75, 0.85. The values of b are extracted from 
the outflow of particles of domains of size up to 250. We 
note that we could obtain an accurate estimate for b in 
this way only for the outflow of the majority species. 
Getting similar estimates for the minority species would 
require significant statistics for much larger domains. 

In Fig. 0|we display the values of 6, as obtained from 
direct measurement of the outflow of majority particles. 
Although a priory one does not expect 6(e, rj) to be given 
by the expression obtained from the ring model of Sec- 
tion HO?] b{e,rj) = 3/2 &r(6, rj), we also display in this 
figure the expression obtained from this formula. We 
find that the numerical data agrees very well with these 
curves. This suggests that in fact the analytical results 
obtained from the ring model (|8T lllf> are valid for non- 
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FIG. 5: Phase diagram for the model, as obtained by setting 
6 = 2 in Eq. |U 

equal densities (77 7^ 1/2) as well. 

We conclude that b{e,rj) = 3/2 6jj(e, 77), where &a(e, 7?) 
is given by ©, provides the correct expression for 6(e, rj) 
for the full model. According to the criterion discussed 
in the introduction, a phase transition into a phase sep- 
arated state is expected at some critical density p for 
b > 2. The transition line in the (e, rj) plane between the 
homogeneous and the phase separated state is depicted 
in Fig. O For values of e which are larger than 0.8 phase 
separation is expected at high densities for any value of 
rj. On the other hand, for e < 0.129 phase separation 
does not occur at any density. 

Direct numerical observation of the predicted transi- 
tion line is hard to obtain. For 77 not too close to 1/2 
one would need to simulate exceedingly large systems, 
far beyond our present reach, in order for the fluid to 
sustain domains which are large enough that the current 
flowing through them takes the asymptotic form. It thus 
remains a challenge to devise some method for numerical 
observation of this transition line. 



IV. SINGLE DOMAIN WITH NON-EQUAL 
DENSITIES: NON-CONSERVING VACANCY 
ENSEMBLE 

As discussed above, for the case of equal densities of 
positive and negative particles one can model a domain of 
length k as an open boundary segment of length k where 
particles enter and exit at the boundaries. This is by 
virtue of the fact that for sufficiently high entry and exit 
rates, the open segment will be maintained in a maximal 
current phase where the bulk density of positive particles 
organises itself to be rybuik = 1/2. One uses the ensemble 
of the open boundary problem to calculate the depen- 
dence of the current on domain length. The maximal 
current phase exhibits long-range correlations that ulti- 
mately generate a slow decay of the current with domain 
length and a condensation transition if e is sufficiently 
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large. 

However an open boundary segment cannot produce 
a bulk density ?7buik 7^ 1/2 and retain long-range corre- 
lations i.e. it cannot produce a maximal current phase 
with X] 7^ 1/2. Thus, for the case of non-equal densities of 
positive and negative particles one cannot use the ensem- 
ble generated by the open boundary problem to calculate 
the current for a domain of size k. Instead one must de- 
vise an alternative ensemble that allows the density of 
the domain to fluctuate about a value r\ 7^ 1/2 whilst re- 
taining the long-range correlations required in a maximal 
current phase. 

In this section we propose an ensemble for the single 
domain of length k that is generated by the dynamics 
of a single vacancy on a ring of size k + 1. The same 
dynamics as the full model is used for the exchange of 
particles (0. Also the dynamics of the vacancy retains 
the processes present in the full model 



+ 
0- 



+ 
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with rate (3 , 
with rate a . 



(15) 



However in addition we introduce two processes where 
particles are not conserved 







0- 



0- 



•0, 



with rate S . 
with rate 7 



(16) 



Note that in (|15|l we have generalised J3J to include a 
rate f3. 

We now show that in the case e = the above dynamics 
generates precisely the ensemble required for a domain 
of length k. To demonstrate this we solve exactly the 
steady state of the non-conserving vacancy problem using 
a matrix product ansatz. 



A. Exact solution for e = 

The matrix product ansatz entails writing the steady 
state as a product of matrices [HI E3 i.e. the steady 
state weight w({si}) for configuration {si} = s%, . . . , Sk- 
is 

w({s i })=Tr[X 1 ---X k ] , (17) 
where the matrix X; is 



X,; = 



Then it can be shown following 0, ^ , that the steady 
state weights for the present model can be written in 
this form provided the matrices D, E and A satisfy the 
quadratic relations 




DE = D + £E 
(3D A = £A , 
aAE = A , 



(18) 
(19) 
(20) 



where £a6 = (3j. For the conserving case 7 = S = 0, £ 
is not fixed and may conveniently be set to £ = 1. This 
recovers the previously known solution 0,E2- However, 
for 7, 8^0, we must take 



£7 
aS 



(21) 



Relations (|18H2(J|) are satisfied if we take A to be the 
projector |V r )(W / |, where we employ a bra-ket notation 
to denote the left and right vectors (W\ and \ V). Then 
letting D = £D relations lfT%H2Uj) reduce to 



DE = D + E , 
0D\V) = \V) , 
a(W\E = (W\ . 



(22) 
(23) 
(24) 



Relations l|22H24|) obeyed by D, E, (W\ and \V), are 
precisely those obeyed by the matrices and vectors used 
to solve the steady state of the open boundary ASEP 
|17| . Thus the weights for the non-conserving vacancy are 
equal to those for an open boundary system reweighted 
by a factor £ N + ; £ acts as a fugacity to tune the relative 
density of positive and negative particles. 

The partition function for the non-conserving vacancy 
system is given by summing over all possible configura- 
tions of positive and negative particles on the ring and 
results in 

Z k = Tr \a(0 + E) k ] = (W\(0 + E) k \V) (25) 

which is precisely the partition function for a domain of 
length k required in the grand canonical partition func- 
tion for the full system (i.e. fixed number of vacancies, 
fluctuating particle numbers, see Section ITD (|4I5I6|I 1 



J M 



-1 M 



Lk=0 



(26) 



Thus we have shown that in the case e = the non- 
conserving vacancy generates the required ensemble for 
domains. For e > we will provide numerical evidence 
that this is still the case. 



B. Numerical Simulations for e > 

We ran numerical simulations of the non-conserving 
vacancy system and measured the density profiles as seen 
from the vacancy. In Fig. 21 we compare the density pro- 
file for the non-conserving vacancy problem at £ = 1 
(corresponding to (r/) = 1/2) with the density profile of 
domains of the same size computed in simulations of the 
full model and those of an open isolated domain. Again 
we find the profiles identical to within the statistical fluc- 
tuations except at sites 1 and k. We have also compared 
the profiles obtained in the non-conserving vacancy en- 
semble at non-equal densities (ry = 0.275) with those ob- 
tained from the direct numerical simulations of the full 
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FIG. 6: Density profile of a domain of size k = 100 at non- 
equal densities rj — 0.275 and e = 1/2. Solid line corre- 
sponds to direct simulation of the full many-domain model. 
Squares correspond to the non-conserving single-vacancy en- 
semble with £ = 1/2. 

model (Fig. [BJ. The very close agreement of the profiles 
in the bulk of the domains provide strong evidence that 
the non-conserving vacancy generates the correct ensem- 
ble for domains. 

To evaluate b(e, rj) we measured the decay of the cur- 
rent Jk(e,rj) with the system size k. As in Section IlIII 
we define the current as the rate at which positive par- 
ticles exchange with the vacancy. The current Jfc(e, rj) is 
obtained by simulating a system of size k + 1, fixing £ 
(Eq. Ell such that the average density of positive par- 
ticles matches rj. The coefficient b(e, rj) is then obtained 
by comparing the currents Jk with the asymptotic form 
^oo (1 + b/k). This requires extensive numerical calcula- 
tion since the current is measured only through the single 
vacancy, taken to be at site 1, and no spatial averaging 
takes place. To overcome this difficulty simulations were 
performed using a multi-spin coding technique p3j . This 
allows many simulations to be run in parallel utilising 
the same random numbers. The resulting b(rf) for vari- 
ous values of e is given in Fig. These results compare 
very well with the analytical results obtained from the 
conserving ring model of Section Til Bl (Eqs. [H HP} ; an d 
with the numerical results of the full model. 



V. GENERALIZATION TO OTHER MODELS 

In previous sections we focused on the model presented 
in Section 2. However, the direct numerical approach for 
calculating b, introduced in Section ITTT1 is applicable for 
other models as well. This method may be used to test 
the assumptions behind the correspondence of a driven 
model to the ZRP, and the applicability of the criterion 
for phase separation, discussed in the introduction. Also, 
in cases where the correspondence between the driven 
model and the ZRP is exact, the question of phase sepa- 
ration in the driven model may be rigorously answered. 

As an example where the numerical method for study- 
ing domains in the full model we consider a two-lane vari- 
ant of the model discussed above. It has been proposed 



that when the rates with e = are considered in 
a two-lane geometry, phase separation may take place 
[24|. In this case the model is defined on a lattice of 
size 2 x L, with periodic boundary conditions in both 
directions, where particles can move cither within a lane 
or between the lanes. Direct numerical simulations of 
the model suggest the existence of phase separation at 
equal densities of the two species. However, calculating 
b by numerical simulations of single open domains yields 
6 ~ 0.8, indicating that phase separation does not occur 
in this model Q . Since no exact mapping to ZRP exists 
in this case, we have tested this result by carrying out 
direct numerical measurement of b in the full model as 
introduced in Section UTTl We find that indeed b ~ 0.8 
also in the full model, verifying the single open domain 
approach in this case. 

Finally, consider another variant of the model of Sec- 
tion 2, whereby the particle exchange rate is replaced 
by H — — > — h with rate 1 and — h — ► H — with rate 
q. This is a generalization of the e = case, allowing 
for backward hopping. This model was introduced in M| 
and studied in 0,SElEi f° r equal densities, and in |25j 
for non-equal densities. With q < 1 this case is qualita- 
tively similar to the e = case considered above, with no 
phase separation taking place at any density. However, 
for q > 1 it can be shown that the steady state weight 
has the same form as (EL with Zk ~ (1 + \^) k q^ k with 
£ = 1 for equal densities 26] . Since this model is exactly 
mapped onto a ZRP, this result may be used to demon- 
strate that the model exhibits phase separation at any 
non-vanishing densities p and rj. From the above result 
for Zk it follows that for large k the current in this case 
is given by Jk = Zk-i/Zk ~ q~ k ^ 2 , which vanishes in 
the limit k — > oo. Thus the criterion correctly predicts 
strong phase separation for all densities of positive and 
negative particles. 



VI. SUMMARY AND DISCUSSION 

In this paper the correspondence between one- 
dimensional two-species driven models and the Zero- 
Range Process is reviewed and extended to consider the 
case of non-equal particle densities. This is demonstrated 
for a two-species exclusion model with 'ferromagnetic' in- 
teractions. To apply this correspondence one has to eval- 
uate the length-dependence of the current emitted from 
a domain of particles. Previous studies were restricted to 
equal densities, were the average velocities of domains is 
zero. In these studies domains were assumed to be statis- 
tically independent, and the current is calculated using a 
model of single domain with open boundaries. In the case 
of non-equal densities domains have a non-zero average 
velocity, and thus this approach is not applicable. 

In the present work we have introduced a method for 
evaluating the current of a domain of length k by direct 
numerical simulation of the full many-domain model. In 
the case of equal densities this method yields the same 



8 



results as before, verifying the validity of the assumptions 
made in formulating the correspondence to ZRP. Namely 
that domains are statistically independent and that the 
current of a domain is given by the steady state current 
of an isolated domain. Moreover, this method may be 
applied to the non-equal density case. 

We also introduced a model for a single domain, which 
enables one to calculate the current of a domain without 
having to resort to a simulation of the full model. Here 
a domain is modelled by a ring with a single vacancy, 
with non-conserving dynamics at the vacant site. It is 
demonstrated that this ensemble yields the same density 
profiles and the currents as domains in the full model. It 
thus provides a rather simple way of analyzing the full 



model. Furthermore, we have outlined an exact solution 
for the non-conserving vacancy model in the case e = 
which extends the range of models solved by the matrix 
product ansatz. 

The phase diagram of the model in the interaction- 
density (e, 77) plane has been calculated using both meth- 
ods, and the transition line to the phase separated state 
has been found. Applications to other models have also 
been discussed. 
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